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SUMMARY OF VIBRATION ANALYSIS 


The objective of this study was to perform a vibration analysis of three 
guyed tower designs for intermediate size wind turbines. These tower designs 
consist of tubular steel sections totaling 28.6 m (94 ft) in length, moxmted 
on a soft pad foundation and supported by guys. Each wind turbine design has 
a 38.1 m (125 ft) diameter rotor and a rated power of 200 to 400 kw with a 

rotor axis at an elevation of 30.5 m (100 ft) above ground. 

This vibration analysis was conducted because there was some concern that 
these guyed towers might have natural frequencies of vibration close to the 
blade passing frequency of the rotor or some multiple thereof. This concern 
arose because of several facts. 

The first is the knowledge that cantilever towers of propeller type wind 
turbines will vibrate violently if excited by the rotor. Secondly, guy cables 
are usually flexible, even with high pretension loads in than. As a result, 

there was concern that the tower frequencies might be too close to the rotor 

speed, Finally it is known that Darrieus type wind turbines, whose central 
column is stabilized by guys attached at the top, have experienced serious 
vibrations. 

Another reason for conducting the vibration analysis was to determine how 
the tower natviral frequencies were affected by the guy dynamic characteristics, 
the location of the guy attachment point on the tower and other factors. 

The method of analysis used is one developed by R. W. Thresher, et al. 
in their report titled "Modeling the Response of Wind Turbines to Atmospheric 
Turbulence", Oregon State University Report No. RLO/2227-81/2 , UC-60, August 
1981. Their report was prepared for the United States Department of Energy, 
Division of Solar Technology, Federal Wind Energy Program. Because of the 
large number of analyses to be performed, a computer program was written in 
FORTRAN which determines the natural frequencies of vibration and plots these 
frequencies as a function of rotor rotational speed. The frequencies of 
vibration determined are for the first two modes of bending-pitch and the 
first two modes of bending-yaw. The prograun determines these frequencies of 
vibration using the inertia and stiffness coefficients of the wind turbine 
system which are calculated using the calculator program. Inertia and Stiff- 
ness Coefficients and Wind Turbines. This calculator program is included 
with the Oregon State University report. 

Two designs in the analysis have two-bladed teetering hub rotors. Since 
the analysis presented by Thresher, et al is limited to three-bladed fixed 
hub rotors, it was necessary to modify their method. The analysis was modified 
by eliminating the rotational inertia of the rotor about its axis and evaluating 
the natural frequencies of vibration at zero rotor speed. For the three-bladed 
fixed hub rotor concept the method was not modified. Furthemore, the method 
of determining the spring rate of the guys used to support the tower had to be 
modified to account for the spring rate of the grouted anchors that the guys 
are attached to. 



In the parametric study several of the system properties were varied to 
determine the sensitivity of the tower's natural frequencies of vibration to 
variances of these properties. Properties varied include the mass of the 
nacelle and rotor, effective guy stiffness, the height of the guy attachment 
point, the location of the nacelle and rotor center of gravity, the mass moment 
of inertias of the nacelle and rotor about the tower axes, and the rotor rota- 
tional speed. Then, by comparing the results of the different concepts, the 
effects of different tower section mass and stiffness properties can be deter- 
mined and the effects of different rotor designs can be evaluated. 

The results showed that only the lowest two frequencies of vibration were 
in an area where they could be excited by the rotation of the rotor. The 
analysis also indicated that these same two modes could be tuned by varying 
the effective guy stiffness, the height of the guy attachment, the mass of the 
nacelle and rotor, or the mass and stiffness properties of the tower sections. 
It was also shown that these same two modes are very close in frequency and 
that both change very little with the rotor's rotational speed. 

In conclusion, the results of the analysis showed that the natural 
frequencies of vibration of all three tower designs are not in a range where 
they can be excited by the rotor. The effects on the natural frequencies of 
vibration caused by changing some of the wind turbine system characteristics 
were also shown. 



INTRODUCTION 


NASA is performing a conceptual design study to investigate the possibility 
of reducing the cost of intermediate size wind turbines. These wind turbines 
are 30.5 m (100 ft) high and have a rated power between 200 and 400 kw. Three 
concepts of tower designs for wind turbines are being considered. Each of these 
tower concepts are made of tubular steel sections mounted on a "soft pad" foun- 
dation and supported by guys with grouted anchors. The "soft pad" foundation 
is a concrete footing with a flexible pad interface between the tower base and 
the footing. This approximates a ball and socket type of joint and minimizes 
bending moments transferred to the footing. Two of the concepts have two-bladed 
teetering hub rotors; the other has a three-bladed fixed hub rotor. All three 
have a rotor diameter of 38.1 m (125 ft) with the axis at 30.5 m (100 ft) 
elevation. 

Cantilever towers of propeller type wind turbines are known to vibrate if 
the tower natural frequencies are close to the blade passing frequency or some 
multiple thereof. Also, in recent years, some towers of Darrieus type wind tur- 
bines which use guys to stabilize the upper end of the tower have experienced 
severe vibrations. These two facts prompted the vibration analysis being re- 
ported here, namely, to determine the vibrational characteristics of the three 
guyed tower concepts under study as mentioned in the previous paragraph. This 
analysis should identify areas where the rotor might excite the tower's natural 
modes of vibration. Furthermore, since these tower designs are still in the 
preliminary design stage, it is important. to determine the sensitivity of the 
natural frequencies of vibration to the various system properties. In this way 
problem areas can be identified and methods of rectification determined. 

A method for vibration analysis of similar wind turbine configurations 
systems has previously been developed by R. W. Thresher, et al of Oregon State 
University. A major part of this analysis is based on their draft report 
titled "Modeling the Response of Wind Turbines to Atmospheric Turbulence", 

Oregon State University Report No. RLO/2227-81/2 , UC-60, August 1981. 

In the analysis of the three tower designs, the various tower and system 
properties of the preliminary designs were estimated and the natural frequen- 
cies determined. The system properties were then varied individually to 
determine their effect on the natural frequencies of vibration. The system 
properties varied include: effective guy stiffness; the height of the guy 

attachment; the mass of the nacelle and rotor; the location of the center of 
gravity of the nacelle and rotor; the mass moment of inertia of the nacelle 
and rotor about the tower's X, Y and Z axes; and the rotor's rotational speed. 
Due to the large number of analyses to be performed, and to expedite the 
calculations, a computer program was written in FORTRAN to calculate and plot 
the natural frequencies of vibration. The results of all the above and the 
computer plots are included in this report. 
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DESIGNS ANALYZED 


Candidate Guyed Tower Designs 

In this report three concepts of tower designs for wind tiarbines were 
analyzed. The three tower concepts have a hi±> height of 30.5 m (100 ft) . One 
of the unique features of these tower designs is that the towers are not can- 
tilevered but are supported by guy rods. These guy rods radiate out from the 
tower at three equally spaced positions. The guys are attached to the tower 
at 10.7 m (35 ft) above the grotind and to concrete pads located 9.1 m (30 ft) 
from the tower's centerline. The concrete pads are anchored by grouted rock 
or soil anchors which extend below the surface to bedrock or to competent soil. 
This is shown in Figure 1. 

To prevent the guys from becoming slack, the guys are pretensioned to a 
value greater than the change in tension experienced under maximum operating 
loads. Similarly, the anchor rods are pretensioned to prevent the anchor pads 
from unseating under all conditions. Thus, the pretension on the anchor rod 
is set to be greater than the maximum tension on the guy rod. 

As mentioned above, the preload prevents the anchor pads from unseating, 
therefore, there is negligible movement of the anchor pads as loads change. 

This causes the effective spring rate of the anchors to be relatively high 
when compared to the spring rate of the guys. 

For the vibration analysis, the effective stiffness of the guys must be 
determined. The derivation of the effective guy stiffness is presented in the 
Thresher, et al report. To include the spring rate of the grouted anchors and 
the soil, the guy system can be modeled as three springs. 
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Thus, the total effective spring rate of the guys and anchors becomes 


2 ^ 
cos 0 


N 


1 

k. 


k + k„ 
a V 


where 

k^ = Total effective guy and anchor spring rate 
k^ = Spring rate of a guy rod 
k^ = Spring rate of an anchor rod 

k^ = Vertical spring rate of soil due to compression 
N = Number of guys per position 
6 = Angle between horizontal guys 

This assumes that the guys and anchors are located at three equally spaced 
positions around the tower and that one guy is attached to one anchor. 

The preliminary design of these towers specifies three guy rods and three 
anchors per position for each concept. Using a guy rod spring rate of 
29.55E + 06 N/m (2.025E + 06 Ib/ft) and a combined anchor and soil spring rate 
of 218. 9E + 06 N/m (15. OE + 06 Ib/ft) , the total effective guy spring rate 
becomes 44.36E + 06 N/m (3.045E + 06 Ib/ft). 

Each tower is made of severoil tubular steel sections; Concept I is made 
of three sections, Concept II has seven, and Concept III has eight. Each of 
these sections has a different wall thickness and, therefore, each have dif- 
ferent inertia and stiffness properties. The lengths and wall thicknesses of 
these sections are listed in Table I and the inertia and stiffness values are 
listed in Table II. The sections are numbered vertically starting at the 
tower base. Additional sections were added to account for the extra mass and 
stiffness of the flanges and the segment from the top of the tower sections to 
the rotor's centerline. 

Located at the top of the tower, as shown in Figures 2 and 3, are the gear 
box, generator, pitch change mechanism and controls. These are all mounted on 
a rotatable platform. The sum of all these parts will be referred to as the 
nacelle throughout this report. Also located at the top of the tower is the 
rotor, which consists of the rotor hub and blades. Concepts I and II both 
have two-bladed rotors with teetering hubs, whereas Concept III has a three- 
bladed rotor with a fixed hub. Furthermore, Concept I has variable pitch 
rotor blades, whereas Concepts II and III have fixed pitch blades. 

The vibration analysis requires the total mass and location of the center 
of gravity for the nacelle and rotor to be calculated. For the various concepts. 
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the total mass and center of gravity locations were estimated as follows : 


CONCEPT 


TOTAL MASS OF 
NACELLE AND ROTOR, 




"'r 


CENTER OF GRAVITY 
LOCATION, q 


I 


16.19E + 03 kg 

(1110 lb secVft) 


2.35 m 
(7.7 ft) 


II 


14.59E + 03 kg 
(1000 lb sec^/ft) 


2.16 m 
(7.1 ft) 


III 


14.18E + 03 kg 
(972 lb sec^/ft) 


-91.43E - 03 m 
(-0.3 ft) 


Since a fixed hub rotor blade can be located much closer to the tower 
than can a teetering hub rotor blade, the rotor of Concept III is located only 
1.52 m (5 ft) from the outside of the tower wall while Concepts I and II have 
their rotor located 3.96 m (13 ft) from the tower. This accounts for the large 
difference in center of gravity locations listed above. The center of gravity 
location is measxared from the tower centerline and a negative value indicates 
that the center of gravity is located downwind of the tower centerline. 


The mass moment of inertia of the nacelle and rotor about the tower axes 
were also calculated, including the mass moment of inertia of the rotor about 
its spin axis. Concepts I and II have teetering hubs; this allows the rotor 
to remain parallel to its plane of rotation during small angular changes in 
pitch and yaw of the nacelle. Thus the rotational inertia about the X and Z 
axes is greatly reduced on Concepts I and II. 

The theory that this analysis is based on does not account for the 
flexibility of the baldes, therefore, it is most applicable to small systems 
where the mass of the blades is a small percentage of the total mass. In this 
analysis the mass of the blades was about 10% of the total mass for Concepts 
I and II, and 13% for Concept III. Also, approximately half the mass of the 
blade was located within the first third of its span. Thus, this simplifica- 
tion is justified. Furthermore, this theory is based on a three-bladed fixed 
hub rotor, therefore, care should be taken in interpreting the results for 
Concepts I and II since they have two-bladed teetering hub rotors. 

The total mass moments of inertia were calculated for each concept and 
are as follows; 


CONCEPT 


MASS MOMENTS OF INERTIA 


kg m/^ 

(lb ft sec^) 


^xx„ 
kg m^ 

(lb ft sec^) 


“■zz. 


kg m 

(lb ft sec^) 


I 


113. 8E + 03 
(83. 9E + 03) 


113. 5E + 03 
(83. 7E + 03) 


169. 8E + 03 
(125. 2E + 03) 
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CONCEPT 


MASS MOMENTS OF INERTIA 


II 


kg m'^ 

(lb ft sec^) 

156. 6E + 03 
(115. 5E + 03) 


^xx- 

kg 

(lb ft sec^) 

125. 4E + 03 
(92. 5E + 03) 


(lb ft sec^) 

132. 3E + 03 
(97. 6E + 03) 


III 


164. lE + 03 
(121. OE + 03) 


133. OE + 03 139. 9E + 03 

(98. IE + 03) (103. 2E + 03) 


Ij. = Rotational mass moment of inertia of the rotor about its 
spin axis (Y) 

Ixx = Mass moment of inertia of the nacelle and rotor about the 
tower ' s X axis 

1^2 = Mass moment of inertia of the nacelle and rotor about the 
tower's Z axis 


The calculator program for determining the inertia and stiffness 
coefficients also requires the number of intervals per tower section to be 
used for integration. The following number of intervals were chosen. 


CONCEPT 

NUMBER OF INTEGRATION 
INTERVALS PER SECTION 

NUMBER OF 
SECTIONS 

TOTAL NUMBER OF 
INTEGRATION INTERVALS 

I 

10 

3 

30 

II 

4 

7 

28 

III 

4 

8 

32 


This completes the required input data for determining the inertia and 
stiffness coefficients for the tower designs. 


Case Studies 

To determine the effects of varying the properties of. the tower, a 
parametric study was performed. In this study each of the tower properties 
were varied, one at a time, to determine their effect on the natural frequency 
of vibration of the tower. First, the effective guy stiffness was varied from 
25% to 200% of the design value. Next, the height of the guy attachment was 
varied from 2/3 to 4/3 of the design height. Similarly, the mass of the 
nacelle and rotor was increased and decreased by 1/3 and the center of gravity 
location was also changed. Finally, the mass moments of inertia about the 
X and Z axes were increased and decreased. The following table lists all of 
the cases examined. Each case was analyzed on each of the three concepts. 
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CASE 


PROPERTY VARIED 


A 25 % of design effective guy stiffness 

B 50% of design effective guy stiffness 

C 100% of design effective guy stiffness 

D 150% of design effective guy stiffness 

E 200% of design effective guy stiffness 

F Height of guy attachment increased by 1/3 

G Height of guy attachment decreased by 1/3 

H Nacelle and rotor mass increased by 1/3 

I Nacelle and rotor mass decreased by 1/3 

J Distance to center of gravity increased by 1/3 


K Distance to center of gravity decreased by 1/3 


L 

Nacelle and rotor 

inertia 

about 

the 

tower ' 

's 

Z 


axis increased by 

1/3 






M 

Nacelle and rotor 

inertia 

about 

the 

tower ' 

's 

Z 


axis decreased by 

1/3 






N 

Nacelle and rotor 

inertia 

about 

the 

tower ' 

's 

X 


axis increased by 

1/3 






0 

Nacelle and rotor 

inertia 

about 

the 

tower ' s 

X 


axis decreased by 

1/3 







P Both of the nacelle and rotor inertias about the 

X and Z axes increased simultaneously by 1/3 

Q Both of the nacelle and rotor inertias about the 

X and Z axes decreased simultaneously by 1/3 

Case C is considered the norm in this analysis. 

The rotor has a design operating rotational speed of 4.19 rad/sec (40 rpm) . 
Since the natural frequency of vibration of a wind turbine system can vary with 
rotor's rotational speed, each concept was analyzed from 0 to 6.3 rad/sec (60 
rpm) at ten equally spaced points. Since Concept I and II have teetering hubs, 
the natural frequencies of vibration will not change with the rotational speed 
of the rotor. Thus, the frequencies of vibration should be evaluated at zero 
rotor speed. On the other hand. Concept III has a fixed hub which causes the 
natural frequencies of vibration to change with rotor speed. Therefore, 

Concept III should be evaluated at operating speed. 
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ANALYTICAL PROCEDURE FOR VIBRATION ANALYSIS 


Inertia and Stiffness Coefficients of Wind Turbine Systems 

To determine the natural frequencies of vibration of wind turbine systems, 
the inertia and stiffness coefficients of the tower must be calculated. The 
determination of the inertia and stiffness coefficients is explained in detail 
in the Thresher, et al report. 

For displacement of the top of the tower in the y direction, the stiffness 
coefficient is (see Figures 4 and 5) 

^ 2 2 
= / EI(z) ('y^(z)) dz + k^’l'J(a) 


where 

V = Translational displacement in y direction 

L = Length of tower 

EI(z) = Tower bending stiffness 

= Displacement function 

= Effective spring constant of guys 

a = Distance from ground to guy connections 

'1' ' , H"' Represent the first and second derivatives with respect 
to Z 

The stiffness coefficient for displacement due to pitching about the X axis is 


2 2 

k = / EI(z) (r'(z)) dz + k '1' (a) 

XX 0 ^ ^ X 


where 


X = Rotational displacement about x axis at the upper end of 
the tower 

For rotational displacement about the tower's Z axis the stiffness coefficient 
is 


k^^ = f GJ(z) ('i':(z))^dz 
9tP n 
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where 


GJ(z) = Tower torsional stiffness 

(|) = Rotational displacement about Z axis at the upper end 

of the tower 

For displacement due to translation in the y direction (V) is influenced 
by the rotational displacement about the x axis (x) < and vice versa. The 
stiffness coefficient relating these two displacements is 

L 

k „ = / EI(z) r'(z)'l'" ( 2 ) dz + k 'l'„(a)’l' (a) 

XV Q V X c V X 

Finally, the stiffness coefficient for translational displacement along the 
X axis is 

k = k - k^ A 
UU W Vx" XX 

where 

U = Translational displacement in the x direction 

The displacement f-unctions (4') are interpolating functions which relate 
the displacements Within the tower to the displacements at the top of the 
tower. For a tower with a ball joint at its base, these displacement f\inctions 
are 

^^( 2 ) = j (z/L) (3 - z/L)^) 

^ ( 2 ) = z/2 (1 - (z/L)^) 

X 

'{'.(z) = z/L 
<p 

'F^(z) = -3z/L^ 

4'"(z) = -3z/L^ 

X 

f'Cz) = 1/L 

(!) 

The equations for a tower with a ball joint at its base were chosen 
because the tower designs utilize bases mounted on flexible pads which reduce 
bending moments transmitted to the foundation. Similar equations for a 
cantilever based tower can be found in the Thresher, et al report. 
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The inertia coefficients of the tower are expressed in a similar fashion. 
They are 


L 

0 



(z) 




where 


m(z) = Mass per lanit length of the tower 


"'r = 


= Mass of nacelle and rotor 


and 


where 


M = / m(z) (z) 4' (z) dz 

VY „ V Y 


M = / m(z) f (z) dz + I 
XX n X XX 


M^^ = /I (z) (z) dz + I 

(})<}) m ({) zz 


I = Mass moment of inertia of nacelle and rotor about the 

XX 

X axis 

I = Mass moment of inertia of nacelle and rotor about the 
zz 

z axis 


I (z) = Mass moment of inertia per unit length of tower about 
the tower centerline 


also 


% = -<"e V ^ 


where 


q = The distance from the tower centerline to the center of gravity 
of the nacelle and rotor 


Finally 


M = M 
UU W 


2{k„ A ) + (k„ A ) (M 

Vx XX vx Vx XX XX 


I ) 

XX 
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The displacement functions are the same as those used in determining 
the stiffness coefficients. These coordinate functions are used in a finite 
element technique and are explained in the Thresher, et al report. 

The stiffness and inertia coefficients as outlined above can be calculated 
using the calculator program Inertia and Stiffness Coefficients of Wind Turbines 
which is included in the Thresher, et al report. The required input data 
include (see Figure 4) : 

1) The type of tower base. 

2) The height of the tower (L) . 

3) The effective spring constant of the guys (k^) . 

4) The height of the guy connection (a) . 

5) The length of the tower sections (L‘) . 

6) The bending stiffness of the tower section (EI(z)). 

7) The torsional stiffness of the tower section (GJ(z)). 

8) The mass per unit length of the tower section (m(z)) . 

9) The mass moment of inertia per unit length about the tower 

centerline (I (z)). 

m 

10) The number of intervals for integration (N) . 

11) The mass moment of inertia of the nacelle and rotor about 

the X and z axes (I and I ) . 

XX zz 

12) The mass of the nacelle and rotor (m^ + m^) . 

13) The distance from the tower centerline to the center of 
gravity of the nacelle and rotor (q) . 

The output from the above program is 


and 


%U' \x' 

’^UU' ^W' ^Vx' ^xx 
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Once the above mass and stiffness coefficients are calculated they can be 
entered into the equations of motion for wind turbine systems. 

The equations of motion for a wind turbine system, excluding aerodynamic 
forces, are 



To obtain the natural frequencies of vibration, these equations of motion 
have the solution 
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Substituting this into the equations of motion yields 


(k - M (jJ ) 
^ UU UU ' 




_ 

0 


% 


o ' 
1 

o 

> 


0 

o 

-e- 


0 

^0 


0 

__ 








(k - M 0) ) 

ux ux 


(il fiw) 
r 


(k -Mu) 
XX XX 


where 

i = /=T 

0 ) = Natural frequency of vibration 

= Mass moment of inertia of the rotor about its spin axis 
= Rotor's rotational speed 

This requires that the determinate of the square matrix be equal to zero. 
Expanding the determinate yields 


UU ffiffl UU <p<j) <p® UU UU d 


" 'Nx’^w 


"vx' - ‘\x"w " '=w“xx " 'VVx' " '"xx^v * \x^> 


(I ^ [(k k ) 

r UU W 
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If we let 


^ ^(j)(j)\u^ 0) + (M 

^ ^XX^VV ^Vx ^^XX^ "'’ ^vv^xx 
^^V^x “ 

- ^’"uu^v + ^w\u 

The expanded determinate becomes 


UU^H “ ^(|) ^ ^ 

"*=vx”vx’ " 


) (m^) + (MyyM^) (O)^)^] 




H „(m^) - 

xv 


“ov'“"> 


0 


Natural Frequencies of Vibration 

At zero rotor speed (Q = 0) , the above becomes 

G,,,(m^) H „{o)^) = 0 

(j)U XV 

thus the roots of and H will give the natural frequencies of vibration 
at taro speed. ♦" *'' 

The equation for can be factored to yield 






2> „ 2, 

W " "W ^ 


Rearranging the determinate 


9 9 

'"ou - V“ "■‘w - «w“ > = ^ 


o2 T 2 2 

U I 0) 

r 


At infinite rotor speed this becomes 

^ V - = ° 

The roots of these will give the natural frequency of vibration at infinite 
rotor speed. 

The calculator program, Natural Frequencies of Wind Turbine Systems, 
included in the Thresher, et al report was written to determine the frequencies 
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of vibration at zero and infinite rotor speed. It can also be used to 
calculate the rotor speed at which a given frequency of vibration will occur. 

In this analysis it was desirable to determine the natural frequencies 
of vibration as a function of rotor speed. Therefore, a computer program was 
written in FORTRAN that extracts the roots in of the expanded determinate. 

=+u'""> ° 


Computer Program for Vibration Analysis 


The FORTRAN program uses the inertia and stiffness coefficients calculated 
on the calculator program, Inertia and Stiffness Coefficients of Wind Turbines, 
and determines the coefficients of ^ and In 

the program 


AG = 
BG = 
CG = 
AH = 
BH = 
CH = 
AD = 
BD = 
CD = 
DD = 


((>()) ^(})c()^UU^ 

2 

M M - M 

UU (})()) U(t) 

2 

^xx^w “ 

Vxx - 
0 

Ww 

-“"h' <%o“w ’'w“uu> 

“uo"w 


Substituting these into the expanded determinate and collecting terms yields 


(AG) (AH) + [(AG)(BH) + (BG) (AH) - 
(CG) (AH) - (CD)](o)^)^ + [(BG)(CH) 
(CG) (CH) (0)^)“^ = 0 


(BD)]0)^ + [(AG)(CH) + (BG) (BH) 
+ (CG) (BH) - (DD) ] (co^) ^ + 


+ 
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After summing like coefficients, the subroutine GNEWTN is called. This 
routine determines the roots of a polynomial using a modified Newton method 
and is capable of locating complex roots. This routine was supplied by NASA - 
Lewis Research Center and requires the subroutine SYNDIV for synthetic division. 
The program is also capable of plotting the natural frequencies of vibration as 
a function of rotor speed and uses plotting routines also supplied by NASA - 
Lewis Research Center, 


The computer program requires the following input data: 

1) Number of data sets (J) . 

2) Number of roots to be plotted (JPO) . 

3) Title. 

4) Mass moment of inertia of the rotor about its spin axis (I^) . 

5) The lowest rotor speed at which roots are to be found (R . ) . 

min 

6) The highest rotor speed at which roots are to be found (R ) , 

max 

7) The number of points between the lowest and highest rotor 
speeds at which roots are to be found (NS) . 


8 ) 

9) 


%<J' ‘SjV “Vv' \x’ “m' "xx- 
''nu' "w S' ’"xx- 


The output includes: 


1) The natural frequencies of vibration at zero rotor speed 

(i.e., the roots of G,„ and H ,). 

(j)U XV 

2) The natural frequencies of vibration at infinite rotor speed 
(i.e., the roots of . 

3) The natural frequencies of vibration at R . , R , and the 
number of points requested between the minimum anl maximum 
rotor speed. 

4) A plot of the natural frequencies of vibration versus rotor 
spin rate. 

A flow chart of this program is included in Figure 6 and a listing can be 
found in Appendix A. 
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RESULTS OF VIBRATION ANALYSIS 


The appropriate tower properties for each case and concept were entered 
into the calculator program. Inertia and Stiffness Coefficients of Wind Tur- 
bines, presented in the Thresher, et al report. The output of the calculator 
program provides the inertia and stiffness coefficients of the wind turbine 
system. The inertia and stiffness coefficients were entered into the FORTRAN 
program. Natural Frequencies of Wind Turbine Systems, which was developed for 
this analysis. The output of the computer program provides natural frequencies 
of vibration as a function of rotor rotational speed. Listed in Tables III and 
IV are the first four modes of vibration for Concepts I and II evaluated at 
zero rotor speed. Table V lists the same modes of vibration for Concept III 
except they are evaluated at a rotor operating speed of 4.19 rad/sec (40 rpm) . 
The results of the different cases are plotted and can be found in Appendix B. 
These four modes of vibration represent the first and second modes of vibration 
in bending coupled with yaw and the first and second mode of vibration in bend- 
ing coupled with pitch. In these tables the first bending-pitch mode is labeled 
(Mode I) , the first bending-yaw mode is labeled o )2 (Mode II) , the second 
bending-yaw mode is labeled W 3 (Mode III) , and the second bending-pitch mode is 
labeled U 4 (Mode IV) . 

Cases A through E for each concept are plotted in Figures 7 and 8 . These 
show the effect of varying the effective guy stiffness on the frequencies of the 
four modes of vibration. Lines labeled IP, 2P and 3P have also been plotted on 
Figure 7 and represent multiples of the rotor's rotational speed. Figures 9 
and 10 show the results of Cases F and G and indicate the effect of changing 
the height of the guy attachment point. The effect of varying the mass of the 
nacelle and rotor. Cases H and L, is shown in Figures 11 and 12 and the effect 
of moving the nacelle and rotor center of gravity is plotted in Figure 13. 
Finally, Cases L through Q, which show the effect of changing the nacelle and 
rotor inertia about the X and Z axes, are plotted in Figures 14, 15 and 16. On 
all of the above plots is the first bending-pitch mode, W 2 is the first 
bending-yaw mode, 0)3 is the second bending-yaw mode, and 104 is the second 
bending-pitch mode. In all cases, up and U 2 were very close and are plotted 
as one line. 

The computer plot of the natural frequencies of vibration versus rotor 
rotational speed for the design case. Case C, can be found in Figures 17, 18 
and 19. The computer plots of the other cases can be foxind in Appendix B. 

The natural frequencies of vibration of U 2 at infinite rotor rotational 
speed are listed in Table VI. Since Concepts I and II have teetering hubs 
which eliminate the gyroscopic effects encountered with a fixed hub, these 
results are not applicable. Even though these results are apparently not of 
any significance now, their importance will be pointed out later. 

The mode shapes of vibration for a typical wind turbine system can be 
seen in Figures 20 and 21. These figures are excerpts from the Thresher, et al 
report and represent a non-guyed wind turbine tower. 
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DISCUSSION OF RESULTS 


The results of the vibration analysis were analyzed to determine the 
important features and to determine the sensitivity of the wind turbine system 
to the variance of individual parameters. Note that in this analysis, in most 
cases, only one parameter at a time was varied, but in actuality when one para- 
meter is changed so are several others. For example, if the mass of the rotor 
changes, the center of gravity of the nacelle and rotor would move and the mass 
moment of inertia about the tower's X and Z axes would also be affected. There- 
fore, several properties of the wind turbine system would change instead of 
just one. 

While analyzing the results of Cases A through E, it was determined that 
all frequencies of vibration are affected by changing the effective guy stiff- 
ness (see Figures 7 and 8 ) . This was expected because changing the effective 
guy stiffness changes the stiffness of the tower. As expected, as the effec- 
tive guy stiffness decreased, so did the natural frequencies of vibration and 
as effective guy stiffness increased, so did the frequencies. Comparing Fig- 
ures 7 and 8 , we find that the lower two modes of vibration are affected the 
most. 


Points of concern are where a multiple of a rotor's rotational speed 
coincides with a natural frequency of vibration. At these points the rotor 
would excite the mode of vibration and cause the tower to oscillate near reso- 
nance. Since Concepts I and II have two-bladed rotors, the first and second 
multiple (IP and 2P) of the rotor's rotational speed are of importance. Since 
Concept III has a three-bladed rotor, the first and third multiples (IP and 
3P) are of interest. When the lowest modes of vibration occur at frequencies 
above the lower rotor orders, the tower is defined as being "stiff". If these 
lower modes of vibration occurred at a frequency below the first rotor order, 
the tower would be considered "soft". In this analysis, the lowest two modes 
of vibration, in most cases, occurred at frequencies well above the first two 
rotor orders for the two-bladed concepts and the first three rotor orders for 
the three-bladed concept. Therefore, all three tower concepts in this analysis 
proved to be "stiff". 

In analyzing the results of Cases A through E (see Appendix B) , it was 
found that the effective guy stiffness had to be decreased considerably before 
any natural frequency of vibration coincided with a lower rotor order. For 
Concepts I and II the effective guy stiffness had to be decreased to 25% of 
the design value before the lowest natural frequency of vibration would be 
excited by the second rotor order. For Concept III the effective guy stiffness 
had to be decreased to 50% of the design value before the lowest natural fre- 
quency could coincide with the third rotor order. This indicates the tower 
designs are very stiff and that the effective guy stiffness could vary signifi- 
cantly before reaching resonance. The high stiffness can be greatly attributed 
to the use of prestressed anchors and anchor pads. Note that the two second 
modes of vibration, CO 3 and 0 J 4 , occur at relatively high frequencies and, 
therefore, should not be of concern. 
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The results of Cases F and G show the effect of changing the height of the 
guy attaching point. Figures 9 and 10 show that as the height of the attaching 
point is raised, the natural frequencies of vibration increase. This implies 
that raising the guy attaching point increases the stiffness of the system and 
thus increases the frequency of vibration. On the other hand, if the height of 
the guy attaching point is lowered, the stiffness of the system decreases as 
well as the natural frequencies of vibration. In both of these cases, the 
first two modes of vibration were affected the most. The first two natural 
frequencies of vibration of Concept III became excited by the third rotor order 
(3P) when the height of the guy attachment was lowered to 2/3 of the design 
value. The first two natural frequencies of vibration for Concepts I and II 
did not coincide with the first two rotor orders in either case. 

In Cases H and I, the mass of the nacelle and rotor was varied. As the 
mass was increased the frequency of vibration of the first two modes decreased 
and as the mass was decreased, frequency increased. Figure 11 shows that, 
within the range of the analysis, none of the concepts have their first two 
modes excited by the rotor orders of concern. Figure 12 shows that the second 
mode of bending-yaw ( 0 ) 3 ) is greatly affected by the mass of the nacelle and 
rotor but the second mode of bending-pitch ( 0 ) 4 ) is not. 

Changing the distance from the tower centerline to the location of the 
nacelle and rotor center of gravity had little effect on Modes I, II and IV, 
but Mode III was greatly affected (see Figure 13) . This plot shows that as 
the distance to the center of gravity is increased, the frequency of vibration 
of the second mode of bending-yaw also increases. The opposite happens when 
the distance is decreased. 

The plots of Cases L through Q show the effect of changing the inertia 
of the nacelle and rotor about the tower's X and Z axes. As the inertia about 
the Z axis was varied, the second mode of bending-yaw ( 0 ) 2 ) was affected the 
most, but as the inertia about the X axis was varied, the second mode of 
bending-pitch ( 0 ) 4 ) was influenced. These effects can be seen in Figures 14 
through 16. The lower modes showed little change and, therefore, were not 
plotted. 

This analysis has shown that only the natural frequencies of the first 
modes of vibration are in a region of concern, i.e. , in a region where they 
might be excited by the rotor. These first modes of vibration, o)]_ and 0 ) 3 , 
were shown to be most greatly affected by changing the effective guy stiffness, 
the height of the guy attachment point, the tower mass and stiffness properties, 
and the mass of the nacelle and rotor. 

Once the basic wind turbine system is designed, the designer cannot 
practically change the mass and stiffness of the tower nor the mass of the 
nacelle and rotor, but the effective guy stiffness can be varied without a 
great deal of redesigning. The effective guy stiffness can be varied by 
changing the size and/or the momber of guy rods. For these tower designs, 
the minimum number and size of the guys were based on the cyclic loading under 
normal operating conditions and the static loads experienced during hurricane 
conditions. Thus, the practical design variables become, increasing the size 
or the number of guys. The effect of the guys can also be changed by varying 
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the height of the guy attachment point, but this would require significant 
redesigning of the tower sections. 

Finally, the results in Table VI, which lists the natural frequencies of 
vibration at infinite rotor rotational speed, show that ^2 changes very little 
with rotor speed, i.e., 012 at zero rotor speed is approximately the same fre- 
quency at infinite rotor speed. It was also shown earlier that at or near zero 
rotor speed, and 0J2 were also approximately of the same frequency. There- 
fore, the frequency of oji and (1)2 at zero rotor speed is approximately the same 
as the frequency of 002 at infinite rotor speed. As discussed in the Analytical 
Procedure, 0^2 infinite rotor speed can be easily found by calculating the 
square root of ^uq/Mpy. Thus, the natural frequencies of vibration of the 
first modes can be easily approximated with one stiffness and one inertia 
coefficient. 
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CONCLUSIONS OF VIBRATION ANALYSIS 


In conclusion, a vibration analysis was performed on three preliminary 
tower designs for intermediate size wind turbines. A parametric study was per- 
formed to determine the sensitivity of the wind turbine system to variances in 
inertia and stiffness properties. The following conclusions were made as a 
result of the vibration analysis. 

1) Under design conditions none of the tower concepts are excited 
by the rotor at or below the operating speed of 40 rpm. 

2) At present, only the lowest two natural frequencies of vibration 
are in frequency range of concern. It was also found in the 
operating range, that these first two modes are very close in 
frequency. Furthermore, the first bending-yaw mode frequency 
( 102 ) changes very little from zero to infinite rotor speed. 

Therefore, the frequency of vibration of most concern can be 
found by evaluating 0)2 at infinite rotor speed. This is simply: 

“2 ' 

where kyy and Myy are the stiffness and mass coefficients for 
translation in the x direction. 

3) The frequencies of vibration of the lowest two modes are 
significantly affected by the effective guy and anchor stiff- 
ness, the height of the guy attachment, and the mass of the 
nacelle and rotor. These modes are also affected by the 
tower section properties, i.e., tower sections' mass and 
stiffness. Thus, the tower may be tuned to avoid excitation 
of the lower modes by changing the effective anchor and guy 
stiffness, the guy attaching height, the mass of the nacelle 
and rotor, or the tower sections' mass and stiffness properties. 

The most practical being the size and the number of the guys. 

4) The center of gravity location of the nacelle and rotor had 
negligible effects on the lower modes of vibration. The same 
is true for the rotational inertia of the nacelle and rotor 
about the tower's X, Y and Z axes. 

5) The frequencies of vibration of the upper two modes were 
affected in all cases, but were of so high a frequency that 
they are of little concern. The effective guy and anchor 
stiffness mostly affected the second bending-pitch mode 
(Mode IV) , whereas the nacelle and rotor mass and the loca- 
tion of the center of gravity affected mostly the second 
bending-yaw mode (Mode III) . The inertia about the Z axis 
had the greatest affect on Mode III, and the inertia about 
the Y axis mostly affected Mode IV. 
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TABLE I 


TOWER SECTIOJ LENGTHS AND THICKNESSES 


TOWER 

CONCEPT 

I 


II 


III 



SECTION LENGTH 

WALL THICKNESS 

SECTION 

m 

mm 

NUMBER 

(ft) 

(in) 

1 

7.62 

7.94 


(25.0) 

(5/16) 

2 

10.67 

9.52 


(35.0) 

(3/8) 

3 

10.36 

7.94 


(34.0) 

(5/16) 

1 

3.66 

7.94 


(12.0) 

(5/16) 

2 

1.83 

9.52 


(6.0) 

(3/8) 

3 

3.05 

12.70 


(10.0) 

(1/2) 

4 

5.49 

15.88 


(18.0) 

(5/8) 

5 

5.79 

12.70 


(19.0) 

(1/2) 

6 

3.66 

9.52 


(12.0) 

(3/8) 

7 

5.18 

7.94 


(17.0) 

(5/16) 

1 

3.66 

14.29 


(12.0) 

(9/16) 

2 

2.44 

11.11 


(8.0) 

(7/16) 

3 

2.44 

14.29 


(8.0) 

(9/16) 

4 

4.57 

19.05 


(15.0) 

(3/4) 

5 

4.57 

15.88 


(15.0) 

(5/8) 

6 

4.88 

12.70 


(16.0) 

(1/2) 

7 

2. 74 

9.52 


(9.0) 

(3/8) 

8 

3.35 

7.94 


(11.0) 

(5/16) 



TABLE II 


TOWER SECTION INERTIA AND STIFFNESS PROPERTIES 


MASS MOMENT 


TOWER 

CONCEPT 

SECTION 

NUMBER 

BENDING 

STIFFNESS 

ElfZ) 

Nm^ 

(lb ft2) 

TORSIONAL 

STIFFNESS 

GJ(Z) 

Nm^ 

(lb ft2) 

MASS PER 
UNIT LENGTH 
M(Z) 
kg/m 

(lb secVft^) 

OF INERTIA PER 
UNIT LENGTH 
Im(Z) 
kgm^ 

(lb ft sec^) 

I 

1 

1.861E+09 

1.489E+09 

279.9 

44.27 



(4.502E+09) 

(3.602E+09) 

(5.845) 

(32 .65) 


2 

2.226E+09 

1 . 780E+09 

335.4 

54.23 



(5.385E+09) 

(4.308E+09) 

(7.006) 

(39.99) 


3 

1.861E+09 

1.489E+09 

279.9 

44.27 



(4.502E+09) 

(3.602E+09) 

(5.845) 

(32.65) 

II 

1 

1.861E+09 

1.489E+09 

280.3 

44.29 



(4.502E+09) 

(3.602E+09) 

(5.855) 

(32.66) 


2 

2.226E+09 

1.780E+09 

336.0 

52.99 



(5.385E+09) 

(4.308E+09) 

(7.018) 

(39.08) 


3 

2.948E+09 

2.358E+09 

447.1 

70.21 



(7.133E+09) 

(5.706E+09) 

(9.337) 

(51.78) 


4 

3.661E+09 

2.929E+09 

557.6 

87.11 



(8.857E+09) 

(7.086E+09) 

(11.645) 

(64.24) 


5 

2.948E+09 

2.358E+09 

447.1 

70.21 



(7.133E+09) 

(5.706E+09) 

(9.337) 

(51.78) 


6 

2.226E+09 

1.780E+09 

336 .0 

52.99 



(5. 385E+09) 

(4.308E+09) 

(7.018) 

(39.08) 


7 

1.861E+09 

1.489E+09 

280.3 

44.29 



(4.502E+09) 

(3.602E+09) 

(5.855) 

(32.66) 

III 

1 

1.861E+09 

1.489E+09 

247.4 

39 ..09 



(4.502E+09) 

(3.602E+09) 

(5.168) 

(28.83) 


2 

2.588E+09 

2.070E+09 

345.7 

54.36 



(6.262E+09) 

(5.009E+09) 

(7.220) 

(40.09) 


3 

3.306E+09 

2.644E+09 

443.5 

69.45 



(7.998E+09) 

(6.398E+09) 

(9.262) 

(51. 22) 


4 

4.364E+09 

3.491E+09 

589.4 

91.69 



(10.56E+09) 

(8.447E+09) 

(12.309) 

(67.62) 


5 

3.661E+09 

2.929E+09 

494.2 

76.93 



(8.857E+09) 

(7.086E+09) 

(10.280) 

(56.73) 


6 

2.948E+09 

2.358E+09 

396.4 

62.21 



(7.133E+09) 

(5.706E+09) 

(8.278) 

(45.88) 


7 

2.226E+09 

1.780E+09 

296.6 

46.78 



(5.385E+09) 

(4.308E+09) 

(6.195) 

(34.50) 


8 

1.86 IE +09 

1.489E+09 

247.4 

39.09 



(4.502E+09) 

(3.602E+09) 

(5.168) 

(28.83) 
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TABLE III 


NATURAL FREQUENCIES 

OF VIBRATION 

OF WIND 

TURBINE TOWER 

DESIGNS 



CONCEPT I 

AT = 0 



CASE 

CONDITION 

FREQUENCY OF VIBRATION, HZ 

(CPM) 



0)1 

“2 

“3 

“4 

A 

0.25 k 

1.26 

1.26 

9.07 

11.35 


T 

(75.9) 

(75.9) 

(544.2) 

(681.0) 

B 

0.50 k„ 

1.72 

1.72 

9.18 

11.77 


T 

(103.2) 

(103.2) 

(550.9) 

(706.4) 

C (NORM) 

1.0 k 

2.26 

2.28 

9.37 

12.60 


1 

(135.5) 

(136.7) 

(562.1) 

(756.0) 

D 

1-5 k 

2.59 

2.62 

9.51 

13.39 


T 

(155.4) 

(157.6) 

(570.9) 

(803.8) 

E 

2.0 k 

2.82 

2.86 

9.63 

14.15 


1 

(169.1) 

(171.9) 

(578.0) 

(849.1) 

F 

1.33 a 

2.83 

2.88 

9.59 

13.05 



(169.8) 

(173.1) 

(575.6) 

(783.2) 

G 

0.66 a 

1.65 

1.66 

9.17 

11.90 



(99.4) 

(99.7) 

(550.5) 

(714.2) 

H 

1.33(m +m ) 

2.02 

2.04 

11.02 

12.52 


K iN 

(121.1) 

(122.8) 

(661.4) 

(751.4) 

I 

0.66 

2.61 

2.61 

8.30 

12.73 



(156.7) 

(156.8) 

(498,3) 

(763.9) 

J 

1.33 q 

2.22 

2.28 

13.52 

12.60 



(133.5) 

(136.7) 

(811.4) 

(756.0) 

K 

0.66 q 

2.28 

2.28 

7.95 

12.60 



(136.7) 

(137.1) 

(477.1) 

(756.0) 

L 

1.33 I „ 

2.26 

2.28 

7.57 

12.60 


ZZ 

(135.4) 

(136.7) 

(454,5) 

(756.0) 

M 

0.66 I 

2.26 

2.28 

13.68 

12.60 


ZZ 

(135.6) 

(136.7) 

(821.1) 

(756.0) 

N 

1.33 I 

2.26 

2.27 

9.37 

12.01 


XX 

(135.5) 

(136.3) 

(562.1) 

(720.7) 

0 

0.66 I 

2.26 

2.28 

9.37 

13.29 


XX 

(135.5) 

(137.2) 

(562.1) 

(797.8) 

P 

1.33(I^„,I ) 

2.26 

2.27 

7.57 

12 .01 


ZZ XX 

(135.4) 

(136.3) 

(454.5) 

(720.7) 

Q 

0.66(I„„,I , ) 

2.26 

2.28 

13.68 

13.29 


ZZ XX 

(135.6) 

(137.2) 

(821.1) 

(797. 8) 
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TABLE IV 


NATURAL FREQUENCIES 

OF VIBRATION 

OF WIND 

TURBINE TOWER 

DESIGNS 



CONCEPT II 

AT = 0 



CASE 

CONDITION 

FREQUENCY OF 

VIBRATION, HZ 

(CPM) 

A 

0.25 k 

T 

“l 

1.27 

(76.2) 

“2 

1.27 

(76.4) 

“3 

12.04 

(722.3) 

^4 

13.04 

(782.3) 

B 

0.50 k 

T 

1.76 

(105.6) 

1.76 

(105,8) 

12.12 

(727.7) 

13.31 

(798.7) 

C (NORM) 

1.0 

2.39 

(143.5) 

2.39 

(143.6) 

12.29 

(737.4) 

13.86 

(831.5) 

D 

1.5 k^ 

2.81 

(168.9) 

2.82 

(169,3) 

12.42 

(745.7) 

14.39 

(863.8) 

E 

2.0 k^ 

3,12 

(187.7) 

3.14 

(188.5) 

12.54 

(752.9) 

14.92 

(895.5) 

F 

1.33 a 

3.07 

(184.0) 

3.08 

(184.8) 

12.47 

(748.7) 

14.12 

(847.3) 

G 

0.66 a 

1.70 

(102.1) 

1.70 

(102.3) 

12,12 

(727.7) 

13.40 

(804.4) 

H 

1 . 33 

2.16 ’ 
(129.9) 

2.17 

(130.4) 

13.90 

(834.2) 

13.65 

(819.1) 

I 

0.66(m^+inj^) 

2.69 

(161.6) 

2.71 

(162.5) 

11.18 

(670.8) 

14.20 

(852.0) 

J 

1.33 q 

2.37 

(142.5) 

2.39 

(143.5) 

15.53 

(931.9) 

13.86 

(831.5) 

K 

0.66 q 

2.39 

(143.5) 

2,40 

(114.3) 

10.89 

(653.3) 

13.86 

(831.5) 

L 

1-33 

ZZ 

2.39 

(143.5) 

2.39 

(143.6) 

10.14 

(608.7) 

13.86 

(831.5) 

M 

0.66 I„„ 
ZZ 

2.39 

(143.5) 

2, 39 
(143.6) 

16.84 

(1010.8) 

13.86 

(831.5) 

N 


2.39 

(143.3) 

2.39 

(143.5) 

12.29 

(737.4) 

13.37 

(802.1) 

0 

hx 

2.39 

(143.5) 

2.39 

(143.8) 

12.29 

(737.4) 

14.41 

(864.6) 

P 


2.39 

(143.3) 

2.39 

(143.5) 

10.14 

(608.7) 

13.36 

(802.1) 

Q 


2.39 

(143.6) 

2.39 

(143.8) 

16.84 

(1010.8) 

14.41 

(864.6) 
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TABLE V 


NATURAL FREQUENCIES 

OF VIBRATION 

OF WIND 

TURBINE TOWER 

DESIGNS 

CASE 

CONCEPT III AT = 4.19 RAD/SEC 

CONDITION FREQUENCY OF VIBRATION, HZ 

(CPM) 



“1 

“2 

“3 

<^4 

A 

0.25 k„ 
T 

1.28 

(77.3) 

1.29 

(77.6) 

10.58 

(635.1) 

13.51 

(811.0) 

B 

0.50 

T 

1.79 

(107.4) 

1.80 

(108.1) 

10.58 

(635.2) 

13.75 

(825.1) 

C (NORM) 

1-0 k„ 
T 

2.44 

(146.8) 

2.47 

(148.2) 

10.58 

(635.3) 

14.22 

(853.3) 

D 

1-5 

T 

2.90 

(174.1) 

2.93 

(176.1) 

10.59 

(635,4) 

14.69 

(881.4) 

E 

2.0 k^ 

3.24 

(194.9) 

3.29 

(197.5) 

10.59 

(635.4) 

15.15 

(909.1) 

F 

1.33 a 

3.16 

(189.6) 

3.20 

(192.0) 

10.58 

(635.3) 

14.44 

(866.5) 

G 

0 . 66 a 

1.73 

(103.9) 

1.74 

(104.6) 

10.58 

(635.2) 

13. 83 
(830.1) 

H 

1.33(m^+rt^) 

2.22 

(133.6) 

2.24 

(134.6) 

10.58 

(635.3) 

13.99 

(839.8) 

I 

0.66(m^+mj^) 

2.74 

(164,9) 

2.78 

(167.0) 

10.58 

(635.2) 

14.59 

(875.6) 

J 

q+0.52m(+l. 7 ft) 

2.44 

(146.8) 

2.47 

(148.2) 

10.65 

(638.9) 

14.22 

(853.3) 

K 

q-0.52m(-l. 7 ft) 

2.44 

(146.8) 

2.47 

(148.1) 

10.71 

(642.9) 

14.22 

(853.4) 

L 

1-33 I22 

2.44 

(146,8) 

2.47 

(148.2) 

9.22 

(553.2) 

14.22 

(853.0) 

M 

0.66 1^2 

2.44 

(146.8) 

2.47 

(148.2) 

12.81 

(768.6) 

14.25 

(855.0) 

N 


2.44 

(146.5) 

2.47 

(148.2) 

10.58 

(635.2) 

13.72 

(823.4) 

0 

^XX 

2.45 

(147.1) 

2.47 

(148.2) 

10.58 

(635.2) 

14.78 

(886.9) 

P 


2.44 

(146.6) 

2.47 

(148.2) 

9.22 

(553.2) 

13.72 

(823.3) 

Q 


2.45 

(146.6) 

2.47 

(148.2) 

12.82 

(767.6) 

14.80 

(886.9) 
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TABLE VI 



NATURAL FREQUENCIES 

OF VIBRATION, HZ (CPM) 



AT INFINITE ROTOR 

ROTATIONAL VEl^OCITY 


CASE 

CONCEPT I 

CONCEPT II 

CONCEPT III 






2 

2 

2 

A 

1.27 

1.28 

1.29 


(76.2) 

(76.8) 

(77.4) 

B 

1.74 

1.77 

1.80 


(104.4) 

(106.2) 

(108.0) 

C 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 

D 

2.67 

2.85 

2.94 


(160.2) 

(171.0) 

(176.4) 

E 

2.92 

3.18 

3.29 


(175.2) 

(190.8) 

(197.4) 

F 

2.93 

3.12 

3.20 


(175.8) 

(187.2) 

(192.0) 

G 

1.67 

1.71 

1.74 


(100.2) 

(102.6) 

(104.4) 

H 

2.07 

2,19 

2.24 


(124.2) 

(131.4) 

(134.4) 

I 

2.66 

2.73 

2. 78 


(159.6) 

(163.8) 

(166.8) 

J 

2.31 

2,41 

2.47 


(138.6) 

(144.6) 

(148.2) 

K 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 

L 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 

M 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 

N 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 

0 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 

P 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 

Q 

2.31 

2.41 

2.47 


(138.6) 

(144.6) 

(148.2) 
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ALL GUYS 

ARE SHOWN 

BUT 

GUY^-.^ 

ONLY ONE 

ANCHOR IS 

SHOWN 

ATTACHMENT 




RING 




/y ' GUY ROD 


/^ANCHOR 
/ PAD 

TOWER 

FOUNDATION 

■ANCHOR 



BED ROCK 


GROUT 

BULB 


WIND TURBINE TOWER GUY AND ANCHOR ROD SYSTEM 
Figure 1 










WIND TURBINE TOWER 

STRUCTURAL AND DIMENSIONAL CHARACTERISTICS 
Figure 4 
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Q = Rotor rotation rate 
(j) = Yaw angle 
X = Pitch angle 
U = Tower top X 
displacement 
V = Tower top Y 
displacement 



WIND DIRECTION 


TOWER COORDINATE SYSTEM 


Figure 5 


























TOWER DESIGNS FOR INTERMEDIATE SIZE WIND TURBINES 


Natural Frequencies of Vibration 
Versus Effective Guy Stiffness 

Concepts I, II and III 
Cases A through E 



, Concept III 
Concept II 

, Concept I 
Wjf Concept II 

Concept III 
<jj^, Concept I 


50% 100% 150% 200% 250% 

Percent Effect Guy Stiffness, %k^ 


The frequencies of vibration of the second two modes as 
the effective guy stiffness is varied from 25% to 200% 
of the design value, = 44.36E + 06 N/m 


Figure 8 







TOWER i;gSIGNS FOR INTERMEDIATE SIZE WIND TURBINES 


Natural Frequencies of Vibration 
Versus Height of Guy Attaching Foint 

Concepts I, II and III 
Cases F and G 



Concept III 
Concept II 


, Concept I 


j_ Uj, Concept II 


- 


-A- 


A 


Uj, Concept III 


Concept I 


_J 1 I I I 

l/3a 2/3a a 4/3a 5/3a 

Guy Attaching Height, a 


The frequency of vibration of the second two inodes as 
the height of the guy attaching point is varied, 
a = 10.7m 


Figure 10 


Natural Frequency of Vibration, HZ 


TOWER DESIGNS FOR INTERMEDIATE SIZE WIND TURBINES 


Natural Frequencies of Vibration 
Versus Nacelle and Rotor Mass 

Concepts I, II and III 
Cases H and I 



Nacelle and Rotor Mass, M = + M^^ 


The frequencies of vibration of the first two modes as 
the Nacelle and Rotor Mass is varied. 

Concept I M = 16.19E + 03 kg 

Concept II M = 14.59E +03 kg 

Concept III M = 14.18E + 03 kg 

Figure 11 
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Natural Frequency of Vibration, HZ 


TOWER DESIGNS FOR INTERMEDIATE SIZE WIND TURBINES 


18.0 


16.0 


14.0 


12.0 


10.0 


8.0 


6.0 1 

0 


Natural Frequencies of Vibration 
Versus Nacelle and Rotor Mass 

Concepts I, II and III 
Cases H and I 



Nacelle and Rotor Mass M = (Mj^ + 

The frequencies of vibration of the second two modes as 
the Nacelle emd Rotor Mass is varied. 

Concept I M = 16.19E + 03 kg 

Concept II M = 14.59E + 03 kg 
Concept III M = 14.18E + 03 kg 


Figure 12 


Natural Frequency of Vibration, HZ 


TOWER DESIGNS FOR INTERMEDIATE SIZE WIND TURBINES 

Natural Frequencies of Vibration 
Versus Center of Gravity Location 


Concepts I, II and III 
Cases J and K 


12.0 


10.0 


8.0 


6.0 L 
0 


/ Concept II 

/ 



_J i J 1 I 

l/3q 2/3q q 4/3q 5/3q 

Distance to Nacelle and Rotor Center of 
Gravity from Tower Center Line, q 


The frequencies of vibration of the second two modes 
as the location of the center of gravity is varied. 
Concept I q = 2 . 35m 
Concept II q = 2 .16m 
Concept III q = -.09m 


Figure 13 



TOWER DESIGNS FOR IWTERMEDIATE SIZE WIMP TURBINES 


Natural Frequencies of Vibration Versus 
Mass Moment of Inertia of the Nacelle and 
Rotor About the Tower's Z Axis 

t 

Cx>ncepts I, II euid III 
Cases L and M 


\ 


\ 



Concept III 
Concept II 


Concept I 


Concept II 


Concept III 


Concept I 


1/3I2Z 2/3IZZ IZZ V3Izz 5/3l2z 


Mass Moment of Inertia of the Nacelle 
and Rotor About the Tower's z Axis, Izz 


The frequencies of vibration of the second two modes 
as the mass moment of inertia of the Nacelle and 
Rotor about the tower's Z axis is varied. 


Figure 14 




TOWER DESIGNS FOR lOTERMEDIATE SIZE WIMP TURBIKES 


Natural Frequencies of Vibration Versus 
Mass Moment of Inertia of the Nacelle and 
Rotor About the Tower's X Axis 

Concepts 1, II and III 
Cases N and O 



, Concept III 
0 )^, Concept II 

Uj/ Concept II 

, Concept I 


-A 

A 

A- 

-o — 

— o — 

— o- - 


lOj, Concept III 
Concept I 


2 

Concept I Ijjjj = 113. 5E + 03 kg m 

Concept II Ijjjj = 125. 4E + 03 kg 

Concept III Ijjjj = 133. OE + 03 kg m 


Mass Moment of Inertia of the Nacelle 
and Rotor About the Tower's X Axis, I^^^j 


The frequencies of vibration of the second two modes 
as the mass moment of inertia of the Nacelle and 
Rotor about the tower's X axis is varied. 


Figure 1 




TOWER DESIGNS FOR INTERMEDIATE SIZE WIND TURBINES 


Natural Frequencies of Vibration Versus 
Mass Moment of Inertia of the Nacelle and 
Rotor About the Tower's 2 and X Axes 

Concepts I, II and III 
Cases P and Q 



1/3 


2/3 


'^zz' '^zz' ’’xx’ 


^^ZZ' ^XX* 


— Concept I 


Uj, Concept II 
ui^r Concept III 


~ (Dj, Concept I 


4/3 


5/3 


*^zz' ^xx^ '^zz' ^xx' 


Mass Moment of Inertia of Nacelle and Rotor 


About the Tower's z and X Axes, d^z' ^xx^ 


The frequencies of vibration of the second two modes 
as the mass moment of inertia of the Nacelle and Rotor 
about the Tower's Z and X Axes are varied simultaneously. 


Figure 16 




imiiiiiimiiiiniiiiiiiiiiiiiiiiiiiiiiiiiiiiiii 

IMIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 

iiiiiiiiimiiiiiiniiiiiiiiniiiiiiiiiiiiimiiii 




“i' “2 


imHBBSII 


l aiiMMm— mjl 

WMMI"' 


4 . 0 O.P. 4 . 5 


R0T0Fi SPEED. RflD/S 

NATURAL FREQUENCIES 0F WIND TURBINE SYSTEMS 

CONCEPT I, CASE C, 100% EFFECTIVE GUY STIFFNESS 

Figure 17 



NfiTURfil FREQUENCY 



NflTURm. FP.EGUENCIES QF rtIND TURBINE STSTFMT. 

CONCEPT II, CASE C, 100% EFFECTIVE GUY STIFFNESS 

Figure 18 






NATURAL FREQUeNCY. HZ 





R8TBR SPEED. RRD/S 

NATURAL FREQUENCIES 0F HIND TURBINE SYSTEMS 


3P 

2 

2P 

IP 


CONCEPT III, CASE C, 100% EFFECTIVE GUY STIFFNESS 

Figure 19 
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BENDING-PITCH VIBRATION MODES AT ZERO ROTOR SPEED 

Figure 20 

(Excerpt from Thresher, et al) 





BENDING- YAW VIBRATION MODES AT ZERO ROTOR SPEED 

Figure 21 

(Excerpt from Thresher, et al) 



APPENDIX A 


FORTRAN COMPUTER PROGRAM 


C*FREE VIBRATION APPROXIMATIONS FOR HORIZONTAL AXIS WIND TURBINES* 

C NATURAL FREQUENCIES OF WIND TURBINE SYSTEMS 
C 

C THIS PROGRAM EVALUATES THE FREE VIBRATION FREQUENCIES FOR A WIND 
C TURBINE SYSTEM. THE REQUIRED INPUTS ARE INERTIA AND STIFFNESS 
C COEFFICIENTS DERIVED FROM CALCULATOR PROGRAM 'INERTIA AND STIFFNESS 
C COEFFICIENTS OF WIND TURBINES'. ALSO REQUIRED ARE THE MASS MOMENT 
C OF INERTIA OF THE ROTOR ABOUT IT'S SPIN AXIS, THE MINIMUM AND 
C MAXIMUM ROTOR SPEEDS AND THE NUMBER OF ROTOR SPEEDS BETWEEN MINIMUM 
C AND MAXIMUM THAT ROOTS ARE TO BE FOUND. 

C 

C THE OUTPUT WILL BE THE ROOTS IN FREQUENCY SQUARED OF THE NATURAL 
C FREQUENCY OF VIBRATION AT ZERO ROTOR SPEED, INFINITE ROTOR SPEED, 

C AND THOSE SPEEDS REQUESTED. 

C 

C THIS PROGRAM WAS WRITTEN FOR NASA LEWIS RESEARCH BY ROBERT J. CHRISTIE 
C OF W. L. TANKSLEY BASED ON WORK OF R.W. THRESHER AND C.E. SMITH OF 
C OREGON STATE UNIVERSITY. 

C 

C TASK ORDER 152-01 FOR NASA CONTRACT NAS3-21900 
C 

0 ************ 

c 

C THIS PROGRAM REQUIRES SUBROUTINE GNEWTN AND SYNDIV 
C 

C R = ROTOR ROTATION SPEED 
C W = NATURAL FREQUENCY OF VIBRATION 
C PLOTTING ARRAYS 

DIMENSION ZZ(200) ,P(400) ,X1(400) 

DIMENSION CX(3) ,CY(4) ,CT{8) 

DIMENSION TITLE (12) 

C 

COMMON /COEFS/ NDGREE,NC,NSF,EPS1,X,COEF(40) ,ROOT(40) 

COMPLEX COEF,X,ROOT 
C PLOTTING DATA 

DATA NX/18/ 

DATA CX/' ROTOR SPEED, RAD/S'/ 

DATA NY/21/ 

DATA CY/' NATURAL FREQUENCY, HZ'/ 

DATA NT/43/ 

DATA CT/' NATURAL FREQUENCIES OF WIND TURBINE SYSTEMS'/ 

C NUMBER OF SIGNIFICANT FIGURES IN ROOTS DESIRED 
NSF=4 

C EPSILON TEST 
EPS1=.001 
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C NUMBER OF SETS OF INPUT DATA AND PLOT OPTION 
READ(5,10) J,JPO 
10 FORMAT (216) 

CALL MODESG(ZZ,12) 

C TITLE 

120 READ(5,20) (TITLE(L) ,L=1,12) 

20 FORMAT (12A6) 

C MOMENT OF INERTIA, MINIMUM AND MAXIMUM ROTOR SPEEDS, NUMBER OF STEPS 
C BETWEEN MINIMUM AND MAXIMUM ROTOR SPEED 
READ(5, 30) RI,RMIN,RMAX,NS 
30 F0RMAT(3E12.6,I6) 

C MASS COEFFICIENTS 

READ(5,40) rmuu,rmuo,rmw,rmvx,rmoo,rmxx 
40 FORMAT(6E12.6) 

C STIFFNESS COEFFICIENTS 

READ(5,50) RKUU,RKW,RKVX,RK00,RKXX 
NP=NS+2 

50 FORMAT(5E12.6) 

C WRITE TITLE 

WRITE(6,180) (TITLE(L) ,L=1,12) 

180 FORMAT (2H1 ,12A6) 

C WRITE MASS MOMENT OF INERTIA AND ROTATIONAL SPEED LIMITS 
WRITE(6,130) RI,RMIN,RMAX 
130 FORMAT (26H MASS MOMENT OF INERTIA = ,E12.6/ 

126H MINIMUM ROTATION SPEED = ,E12.6/ 

126H MAXIMUM ROTATION SPEED = ,E12.6///) 

C WRITE MASS COEFFICIENTS 
WRITE (6, 160) 

WRITE (6, 140) RMUU,RMUO,RMW,RMVX,RMOO,RMXX 
140 FORMAT(7H MUU = ,E12.6/7H MUO = ,E12.6/7H MW = ,E12.6/ 

17H MVX = ,E12.6/7H MOO = ,E12.6/7H MXX = ,E12.6///) 

C WRITE STIFFNESS COEFFICIENTS 
WRITE (6, 170) 

170 FORMAT ( 2 3H STIFFNESS COEFFICIENTS//) 

WRITE (6,150) RKUU,RKW,RKVX,RKOO,RKXX 
150 FORMAT (7H KUU = ,E12.6/7H KW = ,E12,6/7H KVX = ,E12.6/ 

17H KOO = ,E12.6/7H KXX = ,E12.6///) 

160 FORMAT (18H MASS COEFFICIENTS//) 

C ROTOR INERTIA SQUARED 
RISQ=RI**2 

C COEFFICIENTS OF POLYNOMIALS 
C *GOU* 

AG=RKUU*RKOO 

BG= ( RKUU*RMOO+RKOO*RMUU) * ( -1 . ) 

CG= RMUU * RMOO - RMUO * * 2 
C *HXV* 

AH=RKXX*RKW-RKVX* *2 

BH=- ( RKXX*RMW+RKW*RMXX-2 *RKVX*RMVX) 

CH=RMW*RMXX-RMVX* *2 



C FOR ZERO ROTATION SPEED 
COEF(3) =AG 
C0EF(2) =BG 
COEF(l) =CG 

C DEGREE OF POLYNOMIAL 
NDGREE=2 

C NUMBER OF COEFFICIENTS 
NC=3 

C SOLVE FOR ROOTS 
WRITE (6, 190) 

190 FORMAT (20 (2H *) ) 

CALL GNEWTN 
WRITE (6,60) 

60 FORMAT (50H ROOTS OF GOU IN W SQUARED AT ZERO ROTATION SPEED) 
C FOR ZERO ROTATION 
COEF(3) =AH 
C0EF(2)=BH 
COEF(l) =CH 
WRITE(6,190) 

CALL GNEWTN 
WRITE (6, 70) 

70 FORMAT (50H ROOTS OF HXV IN W SQUARED AT ZERO ROTATION SPEED) 
C FOR INFINITE SPEED 
W2=RKUU/RMUU 
W3=RKW/RMW 
WRITE (6, 190) 

WRITE(6,80) W2,W3 

80 FORMAT (34H W SQUARED AT INFINITE ROTOR SPEED/ 2E12.6) 

R=RMAX 

RS= ( RMAX-RMIN) / (NS+1) 

GO TO 100 

110 R=R-RS 
C COEFFICIENTS OF DUV 

100 OMSQ=R**2 

BD=OMSQ* RISQ*RKUU*RKW 
CD=-OMSQ*RISQ* ( RKUU*RMW+RKW*RMUU) 

DD=OMSQ*RISQ*RMUU*RMW 
C COEFFICIENTS OF ENTIRE POLYNOMIAL AT R 
C0EF(5) =AG*AH 
COEF(4) =AG*BH+BG*AH-BD 
COEF ( 3) =AG*CH+BG*BH+CG*AH-CD 
COEF(2) =BG*CH+CG*BH-DD 
COEF(l) =CG*CH. 

NDGREE=4 
NC=5 
X=0 . 

WRITE (6, 190) 

CALL GNEWTN 
WRITE (6, 90) R 

90 FORMAT ( 2 7H ROOTS IN W SQUARED AT R = ,E12.6) 
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C BUILD VECTORS 
M=NS+2 
P(M)=R 

XI ( M) =REAL ( CSQRT ( ROOT (l)))/6.2832 

M=M+NP 

P(M)=R 

X1(M) =REAL( CSQRT (ROOT (2) ))/6.2832 

M=M+NP 

P(M)=R 

Xl(M) =REAL( CSQRT (ROOT (3) ) )/6.2832 

M=M+NP 

P(M) =R 

Xl (M) =REAL( CSQRT (ROOT (4) ) ) /6.2832 
C CHECK NUMBER OF ROTATION FREQUENCY STEPS 
NS=NS-1 

IF (NS. GT. -2) GO TO 110 
WRITE (6, 190) 

WRITE(6,220) 

220 FORMAT ( 5 5H ROTOR, RAD/S Wl HZ W2 HZ W3 HZ W4 HZ) 

WRITE(6,210) (P(L) ,Xl(L) ,Xl(L+NP) ,Xl(L+2*NP) ,Xl(L+3*NP) ,L=1,NP) 

210 F0RMAT(5(E12.6) ) 

IF( JPO.EQ.O) GO TO 230 

IF(JPO.EQ.l)N=NP 

IF(JPO.EQ.2)N=2*NP 

IF(JPO.EQ.3)N=3*NP 

IF(JPO.EQ.4)N=4*NP 

N=N+1 

Xl(N)=0. 

P(N)=0. 

CALL GRAPHG ( ZZ ,N , P , Xl , NX , CX ,NY , CY ,NT , CT) 

CALL PAGEG(ZZ, 0,0,1) 

230 J=J-1 

IF(J.GT.O) GO TO 120 
CALL EXITG(ZZ) 

STOP 

SUBROUTINE GNEWTN 
C 

C POLYNOMIAL ROOT FINDER, USES A MODIFIED NEWTON METHOD 
C INPUT TO THIS SUBROUTINE IS AS FOLLOWS 
C NDGREE DEGREE OF POLYNOMIAL INTEGER 

C NSF NO. OF SIGNIFICANT FIGURES ACCURACY DESIRED (INTEGER) 

C EPSl AN EPSILON TEST, USUALLY APPROX .0001 

C COEF(I) ,I=1,NC COEFFICIENTS (COMPLEX) 

C CALLING PROGRAM MUST HAVE THE FOLLOWING CARDS 
C COMMON /COEFFS/ NDGREE, NC, NSF, EPSl ,X,COEF(40) ,ROOT(40) 

C COMPLEX COEF,X,ROOT 

C X INITIAL GUESS OF A ROOT (SUB USES X=0 IF NOT SPECIFIED) 

C 
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COMMON /COEFS/ NDGREE,NC,NSF,EPS1,X,COEF(40) ,RDOT(40) 

COMMON /STORE/ C(40) 

COMPLEX COEFFS 
COMPLEX PVAL 

COMPLEX COEF,C,F, DELTA, FPR,X,XO, ROOT 
COMPLEX FO 

DIMENSION PVAL(20) ,NP(21) ,COEFPS(40) 

EQUIVALENCE (AVX,NAVX) , (AVXO,NAVXO) 

LOGICAL SKIPST 

DATA MASK7/0777777777777/ 

N=NDGREE 

NC=N+1 

IF(EPS1 .LE. 0.) EPS1=0.0001 
2 NNSF=(9-NSF) *3 
MASK=0 
NZ=36-NNSF 

FLD ( 0 ,NZ , MASK) =FLD ( 0 , NZ , MASK7) 

MTEST=MASK 
DO 79 1=1, NC 
NP(I) =NC-I 
79 COEFFS(I) =COEF(I) 

WRITE (6, 200) N, (NP(I) ,COEFFS{I) ,I=1,NC) 

200 FORMAT (22H POLYNOMIAL OF DEGREE 13/1 7H WITH COEFICIENTS/ 

1 (5H X** I2,1P2E18.7) ) 

WRITE (6, 204) X 

204 FORMAT (31HK USING STARTING VALUE FOR X= 2F15.8// 

1 4X,1HM,4X,1HK,3X,2HIT,6X,4HABSF,14X,6HABSFPR,12X,1HX,32X,2HX0) 
NR=0 
K=0 
IT=0 
MN=0 

19 CALL SYNDIV(N,X) 

SKIPST=. FALSE. 

F=C(NC) 

AVF=CABS(F) 

IF(AVF .LT. l.OE-10 ) GO TO 50 
M=1 

DO 10 I=N, 32767, 32767 
FPR=C(I) 

AVFPR=CABS(FPR) 

IF(AVF .EQ. 0.) GO TO 50 

IF(AVFPR .EQ. 0.) GO TO 10 

DELTA=-F/FPR 

IF(M .EQ. 1) GO TO 9 

IF ( CABS ( DELTA) .LT. l.OE-14) GO TO 9 

SM=M 

DELTA=CLOG ( DELTA) /SM 
DELTA=CEXP ( DELTA) 

9 IF (CABS (DELTA) .LT. 2 . *AVF) GO TO 11 

10 M=M+1 

11 IF (IT .EQ. 0) GO TO 6 
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IF(AVF .LT. STORE) GO TO 666 
IF(M .NE. MS) GO TO 66 
SKIPST= .TRUE. 

GO TO 16 

666 IF(REAL(F*FO) .GE. 0.) GO TO 6 
SKIPST= .TRUE. 

GO TO 16 
6 XO=X 

AVXO=CABS(XO) 

ST0RE=AVF 

IT=IT+1 

IF(MOD(IT,15) .EQ. 0) MrEST=MTEST*2 
SKIPST= .FALSE. 

K=0 

66 IF(M .EQ. 1) GO TO 15 
H=2**K 

DELTA=DELTA/H 

15 X= X+2.* DELTA 

16 X=(X +XO)/2. 

K=K+1 

WRITE(6,201) M,K,IT,AVF,AVFPR,X,X0 

201 F0RMAT(3I5,6G18.8) 

IF(SKIPST .AND. K .LT. 10) GO TO 19 
AVX=CABS (X) 

N XMXO=N AVX-NAVXO 

IF ( CABS ( X-XO) .LT. l.OE-6 ) GO. TO 50 
IF(AND(NXMXO,MTEST) .EQ. 0.) GO TO 50 
MS=M 
GO TO 19 

50 NR=NR+1 

WRITE (6, 505) NR,X 

505 FORMATdOHL ROOT NO. 14, 9H FOUND, = 2G15.7) 

ROOT (NR) =X 

IF (NR .EQ. NDGREE) GO TO 100 
MINUSN=-N 

CALL SYNDIV(MINUSN,X) 

N=N-1 
NC=NC-1 
. DO 51 1=1, NC 

51 COEF(I)=C(I) 

IT=0 

MTEST=MASK 

X=CONJG(X) 

GO TO 19 
100 CONTINUE 

NC=NDGREE+1 

WRITE (6, 200) NDGREE, (NP( I) ,COEFFS(I) ,I=1,NC) 

WRITE (6, 202) 

202 FORMAT(1HL,20X,32HTHE ROOTS OF THIS POLYNOMIAL ARE / 

1 4HLNO. ,8X,4HREAL,10X,9HIMAGINARY,10X,9HMAGN1TUDE,9X,5HANGLE) 
PROD=l . 



DO 101 I=1,NDGREE 
PV7Ui(I)=COEFFS(l) 

DO 102 J=2,NC 

102 PVAL(I) =PVAL(I) *POOT(I) +COEFFS(J) 

XMAG=CABS ( ROOT ( I) ) 

PROD=PROD*XMAG 

XANG=57.29578*AIMAG(CLOG(ROOT(I) ) ) 

101 WRITE (6, 203) I, ROOT (I) ,XMAG,XANG,PVAL(I) 

WRITE (6,205) PROD 

203 FORMAT(I4,lP3G17.7,OPFll.l,5H DEG. ,1P2G17, 7) 

205 FORMAT ( 2 7HL THE PRODUCT OF THE ROOTS = G20.8) 

9999 CONTINUE 
RETURN 

SUBROUTINE SYNDIV (NDEG , D) 

C 

C SYNTHETIC DIVISION SUBROUTINE 

C N= DEGREE OF POLY, AND IF N.L.T. 0, DO ONLY 1 LINE OF DIVISION 
C D= DIVISOR 

C C= ARRAY OF (N+1) COEFFICIENTS ClX**N+C2X** (N-1) CN 

C 

COMPLEX D,C,COEF,X,ROOT 

COMMON /COEFS/ NDGREE ,NC ,NSF,EPS1 ,X,COEF (40) ,ROOT(40) 

COMMON /STORE/ C(40) 

L=NDEG 

NN=IABS(L) 

IF(L .LE. 0) L=1 
DO 2 1=1 ,NC 
2 C(I)=COEF(I) 

DO 1 1=1, L 
M=NN+2-I 
DO 1 J=2,M 

1 C(J) =C(J-1) *D +C(J) 

RETURN 

END 
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APPENDIX C 


SYMBOLS 



Distance from ground to guy connections 

Cross sectional area 

Modulus of elasticity 

Flexibility influence coefficients 

Loads in u, v directions 

Shear modulus of elasticity 

Hub radius 

Area moment of inertia 

Mass moment of inertia per unit length about the centerline of tower 

Mass moment of inertia about the rotor spin axis 

Mass moment of inertia of nacelle-rotor around x, z axes 

Polar moment of inertia 

Anchor stiffness 

Effective spring constant of guys 

Spring rate of a guy rod 

Effective guy and anchor stiffness 

Vertical spring rate of soil due to compression 

Stiffness influence coefficients 

Guy cable length 

Beam, tower length 

Tower section length 

Mass per unit length 

Tower inertia coefficients 
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m. . 
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O.P. 


P. 
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U, V 


w 
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0 

A 

<!>/ X 
1 

0 ) . 

X 

n 
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Mass of nacelle 
Mass of rotor 

System inertia coefficients 
Total mass of tower 

Number of integration intervals, also number of guys per position 
Operating point 

Generalized load at ith location 
Rotor's rotational velocity 
Distance to nacelle-rotor C.G. 

Thermal energy 
Strain energy 
Kinetic energy 

Displacements at tower upper end 
Work 

Approximate function for blade model deflections 

Angle between guys and ground 

Dimensionless blade natxiral frequencies 

Rotational displacements at tower upper end 

Displacement functions 

Natural vibration frequencies 

Rotor speed 

Dimensionless rotor speed frequency 
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